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TECHNICAL NOTE 3719 

APPLICATION CF SCATTERING THEORY TO THE MEASUREMENT OF 
TURBULENT DENSITY FLUCTUATIONS BY AN OPTICAL METHOD 
By Howard A. Stine and Warren Winovieh 


SUMMARY 


An analysis shows that the scattering of unpolarized plane light 
waves which penetrate turbulent , transparent gases provides a measure of 
the average integral scale and intensity of the turbulent density fluctu- 
ations. The analysis is based upon the scattering coefficient deduced 
by Booker and Gordon and, thus, assumes isotropic fluctuations having an 
exponential] y decaying (Markoff ian) correlation function. Photometric 
data taken through turbulent boundary layers in air are found to conform 
functional] y with the analytical prediction. 


INTRCDUCTION 


The effect of random variations in index of refraction of the atmos- 
phere on the propagation of electromagnetic waves has recently received 
widespread theoretical consideration (refs. 1 to 3 ) . Attention, in the 
main, has been concentrated on the problem of the "forward scattering" 
of vhf and uhf radio waves with the intent to explain unexpectedly high 
signal strengths which exist beyond the line of sight. In general, it 
has been demonstrated that this scattering of beamed radiation out of the 
line of sight can be explained by postulating reasonable values for the 
intensity and scale of turbulent density fluctuations in the atmosphere. 

An analogous scattering occurs when a wave of visible light traverses 
a turbulent boundary layer, wake, or jet. The significant interest in the 
scatter phenomenon for light transmission is not reception beyond the line 
of sight, but the examination of the scattered light field with the intent 
of deducing the turbulent parameters that give rise to scattering. 

Liepmann (ref. k) studied the fluctuations in the deflection angle 
of a light ray passing through a turbulent boundary layer and obtained 
the mean square deflection in terms of the correlation function and the 
mean square fluctuations in refractive index. The results at radio wave 
lengths and Liepmann* s analysis suggest that measurements in the scattered 
field can be applied to determine the turbulent properties of the medium. 

The purpose of this report is to describe an optical method for 
measuring the average intensity and integral scale of turbulent density 
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fluctuations in transparent gases. The method relies upon the hypothesis 
that the turbulent field can he uniquely defined in terms of the scattered 
electromagnetic field which arises when a plane wave passes through the 
turbulence. Photometric exploration of the scattered field can then be 
interpreted with the aid of existing theory to yield the turbulence struc- 
ture as characterized by the integral scale and average intensity of the 
density fluctuations. The application of forward scattering theory to 
such measurements was discussed in a paper by Stine and Winovich given 
before the joint meeting of the American and Mexican Physical Societies 
in Mexico City, August 31^ 1955* This report is an expanded version of 
the paper. 

The significance of the theory developed here is that it provides a 
means for measuring turbulent density fields by relatively simple optical 
methods. Although the technique has limi tations, it has the singular 
advantage of responding only to density fluctuations as opposed to the 
mixed-response signal of the hot-wire anemometer. Thus, it promises to 
be especially applicable to turbulence measurements in compressible flows. 


NOTATION 
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velocity of light in a vacuum, ft/sec 
Gladstone-Dale constant, 0.117 ft 3 /slug (eq. (a 6)) 
focal-plane aperture diameter, ft (fig. 2) 
average intensity of light, watts/ft 2 

intensity or transmittance ratio relative to incident wave 
focal length of telephotometer, ft (fig. 2) 
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— J dimensionless grouping indicative of the principal scat- 
tering angle for the Markoff i an correlation function 
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dimensionless grouping indicative of the principal scat- 


tering angle for the Gaussian correlation function 
integral scale of density fluctuations, ft 
distance through turbulence along path of primary beam, ft 
velocity of light in air, ft/sec 
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power associated with electromagnetic wave, watts 


correlation function for turbulent fluctuations, dimensionless 


attenuation coefficient, ft" 1 

absolute dielectric constant, (emu) 2 lb" 1 ft" 4 sec 2 


index of refraction 



scattering angle, radians 

a 

half-angle subtended by focal-plane aperture, tan -1 — 

2F 


s 


radians 


wave length of radiation, ft 
absolute permeability, lb ft 2 (emu) 
density, slug/ft 3 

scattering cross section (steradian ft) 1 , equation (6) 
solid angle corresponding to 8, steradians, ^ = sin 2 ^ 


angle between incident electric vector and the scattered ray 


Subscripts and Superscripts 
Gaussian form of correlation function 

value measured by telephotometer with focal-plane aperture 0 M 

scattered 

total 

Markoffian form of correlation function 
value associated with incident wave 
value pertaining to measurements in a vacuum 
mean value 
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THEORY - 


Consideration is directed first toward an application of the pertinent 
scattering theory to the attenuation of a plane light wave traversing a 
scattering medium. To this end, consider the telephotometer system shown 
schematically in figure 1. A uniform, parallel "beam of monochromatic 
light is directed through a turbulent medium, undergoes scattering, and 
is received by a telephotometer. Tbe scattered field is e x a m i n ed at the 
focus of* the telephotcmeter by means of calibrated apertures used in con- 
junction with a photomultiplier tube. The action of the axially symmetric 
telephotometer of figure 1 is shown schematically in figure 2. The sketch 
illustrates the propagation path for one of the scattered rays that just 
impinges on the periphery of the aperture. By the geometry of the system, 
fll 1 rays scattered through angles greater than 0 y are rejected by the 
instrument. So long as the angle 0 M is sufficiently small, all other 
rays making this angle with the optical axis will also impinge on the 
periphery of the aperture. Thus, it is clear that the angle 0 in this 
case is just equal to the an gnl ar aperture of the telephotometer, 0 ^. 

The instrument with an aperture 0j^, then, records only those rays scat- 
tered through angles that do not exceed 0j^. As the aperture size 
increases, the telephotometer receives a greater percentage of the scat- 
tered power, and the transmittance ratio derived from observations with 
the telephotometer increases. 

For a given focal-plane aperture, the attenuation that occurs as the 
light beam traverses the turbulent medium is governed by the considerations 
which follow. It is known that aggregations of particles such as air, 
smoke in air, and aerosols attenuate a beam of parallel monochromatic 
radiation in close accord with Lambert’s exponential law for homogeneous 
media (refs. 5 , 6 , and 7 ), which states that the percentage of intensity 
loss is proportional to the distance traversed: 

§ —a ay ( 1 ) 

Tht* attenuation per foot of penetration, a, is dependent on the wave- 
length and state of polarization of the incident radiation, the structure 
of the scattering particles, and the intensity of fluctuations in particle 
number density. Ordinarily, a is considered as consisting of two dis- 
tinct parts, an absorption coefficient and a scattering coefficient, to 
distinguish between the intensity decrease due to disappearance of visible 
light into radiation at other frequencies (such as heat) and that due to 
deflection of visible light from the direction of primary propagation. 
Because visible light is absorbed only very slightly by air, the attenua- 
tion per foot in equation (l) can be taken as that due to scattering of 
the radiation. 

An Important consequence of the exponential form of Lambert’s law 
is that measurements madg on two or more identical scattering layers 
can be reduced to a unit penetration. If the radiation passes through 
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n identical layers (or is reflected "back along the same path n times), 
the logarithmic solution of equation (l) for a constant mean scattering 
coefficient relates the energy loss for a unit penetration to the total 
loss for n units of penetration by the expression 



where e/E q is the transmittance ratio and the subscripts refer to the 
ratio of the path lengths. For example, if the radiation passes through 
two identical layers and is measured, the equivalent loss for one layer 
is obtained by taking the square root of the measured quantity. 

Consideration is directed, secondly, to the effects of the turbulent 
density fluctuations on the attenuation of a light beam that obeys Lambert* 
law. The structure of the turbulence affects the transmittance ratio by 
entering into the specification of the attenuation or scattering coeffi- 
cient, a. Light, considered as an electromagnetic wave, consists of an 
electric field and a magnetic field which are mutually perpendicular to 
each other and to the direction of propagation. As the light wave trav- 
erses a medium in which the particle number density fluctuates, the inter- 
action of the electrons associated with the particles and the electric 
field of the wave causes a reradiation of energy to occur. The reradia- 
tion is a manifestation of the accelerations of the light wave. The 
magnitude of the intensity of the scattered rays is expressed as a frac- 
tional part of the primary beam intensity 

P s = *Eo (3) 

The ratio of the scattered power to the incident intensity is denoted by 
a, termed the classical scattering coefficient. By definition of ct, the 
incident radiation falling on a unit area of a surface perpendicular to 
the primary beam loses a fraction of its energy equal to a. Thus, in 
unit time the medium scatters as much radiation as falls on an area equal 
to a; and, for this reason, o is commonly referred to as the classical 
cross section for scattering. 

In passing through a scattering volume, an electromagnetic wave 
undergoes scattering in all directions. This physical result is embodied 
in the concept of the scattering cross section described previously. 
Specifically, the scattering cross section is defined as the fraction of 
the incident radiant flux scattered per unit solid angle per unit volume. 
When a volume containing a turbulent medium is viewed over the solid angle 
ux to ua, the average attenuation coefficient is found by integrating the 
scattering cross section over the solid angle range: 
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where 0 ( 10 ) represents the scattering cross section in terms of the solid 
angle w. For line-of-sight propagation and axially symmetric receiving 
systems , the solid angle to is related to the scattering angle 0, meas- 
ured from the direction of primary propagation, "by the identity 

■— s sin 2 4 (5) 


The scattering cross section can he found "by solving Maxwell’s wave 
equations in a region characterized by a fluctuating dielectric constant 
e = e + As with the Ae’s correlated in space by a correlation function 
B(y/l). Booker and Gordon (ref. 1) obtained a specification of the scat- 
tering cross section for a turbulence model with isotropic fluctuations 
and a space correlation that obeys an exponential (Markoffian) decay law. 
Their result is: 


o(e) » 


(“T w 


sin^X 


For forward scattering 1 - which occurs at nearly right angles to the 
incident electric field - the term sin 2 X can be taken as unity. This 
simplifies the scattering cross section so that 




From equation (7) , it can be seen that the power distribution in the 
scattered field is symmetrical about the propagation axis and is propor- 
tional to the mean square of the dielectric fluctuations, (Ae/e) 2 , and 
is also directly proportional to the cube of the integral scale, Z. 
Equation (7) indicates that the transmittance ratio will be considerably 
affected by relatively small changes in the turbulence structure. The 

1 Rigorously, the forward scattering assumption need not be intro- 
duced here. Equation (6) is in a form applicable to polarized radiation; 
but it can be modified to be applicable to unpolarized radiation by 
replacing the term sin 2 X by 1 + cos 2 0/2 (ref. 8). For the small angles 
associated with forward scattering, the latter term is also essentially 
unity. In the subsequent development, the results are unaltered by the 
simplification employed here. 
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turbulent density fluctuations are simply related to fluctuations in the 
dielectric constant as shown in Appendix A. 

The simplification of equation (6) to equation ( 7 ) by restricting 
the discussion to forward scattering can be justified by considering the 
effects on the cross section of the ratio bit l/~K. By integration of equa- 
tion (7) over the limi ts zero to an arbitrary scatter ing angle 0, the 
scattered power found within the cone generated by the corresponding Bolid 
angle can be obtained. The result of such an integration is: 

1 - sln2 | 

1 + 1 

The half-power cone, then, occurs for the angle that has 




In anticipation of the result, the sine can be replaced by the angle itself 
and the factor 2 neglected to define the principal scattering cone as one 
with a half- angle given by: 


9 


A 

2x1 


Upon taking l/lO inch as an order of magnitude of the integral scale in 
a turbulent boundary layer and using 2x10" 5 inch (5200 A) for the wave 
length of light, one finds that the half-power cone has a half- angle of 
about 2x10" 5 radians (about 4 seconds of arc). This justifies the assump- 
tion that the scattered wave is beamed mainly in the forward direction. 

The power polar diagram of the scattering cross-section function shown in 
figure 3 represents the scattering from a unit volume of turbulence and 
illustrates the corresponding half-power cone described previously. 


A general expression can now be derived relating the ratio of the 
measured intensities that result when two arbitrarily sized apertures 
that subtend solid angles and 0 ^ (wa > o^) are placed in the focal 
plane of the telephotometer of figure 1. Equation (1) can be integrated 
in closed form by utilizing the scattering cross-section expression given 
by equation (k ) . Integrating (l) between the limits (jy to Ug and then 
0 to y gives 



a dw dy 


(8) 
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When "the second aperture, cjg , is selected to "be much larger than the aper- 
ture that receives the half -power cone, virtually the entire scattered 
field can he received. For this condition, the upper limit (in 
eq. (8)) can be relaxed to kit; and E s can be replaced by E 0 , the inten- 
sity associated with the incident wave. In practice under conditions of 
forward scatter ing , this reference aperture is act ually small because the 
half -power cone itself is small. The integration can be readily accom- 
plished with the aid of the identity relating and (eq. (5)) and 
the scattering cross section of equation (7). After integrating and using 
the assumption that the integral scale of the fluctuations is independent 
of depth of penetration, one obtains the equation for the telephotometer: 


, % 

In — 

E n 


dr- 


sin 2 


(JL + K sin 2 


where 


K 


<¥J 


With some rearrangement, this becomes 


(dhc * 0‘ 



(9) 



This form of the telephotometer equation relates the transmittance ratio 
E m /E 0 measured by the telephotcmeter with a focal-plane aperture 9^ to 
the turbulence structure characterized by the integral scale and average 
intensity of turbulent fluctuations. Substantially the same result occurs 
for the assumption of a Gaussian form for the correlation function 
(Appendix B) . 

The form of the telephotcmeter expression (eq. (10)) is significant. 
As shown in figure 4, a plot of the transmittance ratio function 
-l/Zn(Ejij/E 0 ) is linear with the aperture function tan 2 (e M /2). On the 
linear plot, the mean structure of the turbulence is contained in the 
slope and intercept - or graphically in the x and y axis intercepts. 

The x and y axis intercepts define the integral scale and average 
intensity of turbulent fluctuations, respectively. Thus, a series of 
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attenuation measurements obtained from photometric surveys of the light 
field of a given wave length through apertures of varying diameter can 
yield values of the unknown turbulent parameters. 


EXPERIMENTAL RESULTS 


Figure 5 shows the extent to which the foregoing analysis is borne 
out by experiment. The telephotometer arrangement was similar to that 
shown schematically in figure 1; and the density fluctuations occurred 
in the two turbulent boundary layers on the side walls of the 1- by 3~foot 
supersonic wind tunnel no. 1 at the Ames Aeronautical Laboratory. Free- 
stream conditions were maintained constant and the boundary layers were 
altered from the normal smooth-wall case both by means of single roughness 
elements located upstream of the throat and bands of distributed roughness 
of different extent along the tunnel side walls. Transmittance ratios 
measured with a set of focal-plane apertures are plotted in the functional 
form suggested by equation (10) for each boundary- layer configuration. 

The linear result predicted by the analysis was obtained. For the condi- 
tions shown, the integral scales range from 1/10 to l /25 of the boundary- 
layer thickness, and the average fluctuation intensities range up to about 
5 percent of the mean boundary- layer density. Analogous values of 2/5 and 
5 percent for velocity turbulence have been measured in low- speed flows by 
Klebanoff and Diehl (ref. 9) • 


DISCUSSION 


Over the range of flow conditions tested, the telephotometer data 
exhibit the functional dependence between transmittance ratio and aperture 
function predicted by the analysis. This experimental agreement lends 
weight to the assumptions made in the analysis and allows interpretation 
of the x and y axis intercepts of equation (10) (fig. *0 in terms of 
the turbulence parameters. However, the interpretation placed upon the 
intercepts is influenced by the turbulence model used in the analysis. 

In particular, the interpretation is affected by the mathematical form 
of the correlation function. 


The specification of the scattering cross section (eq. (6)) is 
directly influenced by the form of the correlation function. The scat- 
tering cross section used in this paper is based upon the Markoffian 


correlation function 




This form was selected because 


it has been found applicable to forward scattering phenomena in the radio- 
frequency portion of the spectrum. It should be noted that this form is 
not necessarily the best for application to light scattering by turbulent 
boundary layers. 
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Physically, the correlation function obtained by measurements is 
characterized by zero slope at the origin which can be ascribed to vis- 
cous damping over a region characterized by the microscale of turbulence 
(ref. 10). The Markoff ian function does not possess this property. A 

Cauchy form (z(y/l) = l/[l+ (y/l) 2 ] 2 ) or a Gaussian form ^R(y/l) = e"^ y / 1 ^ ^ 

of correlation function is physically better suited because the slope is 
zero at the origin for these functions. Using the cross section expres- 
sions obtained unde r these two assumptions (Wheelon and Muchmore , ref. 3 , 
p. 1 ^ 52 ), one can obtain the equation for the telephotometer for these two 
forms. For the Cauchy form, the result leads to a functional dependence 
tha t is not verified by the optical data. For the Gaussian form, however, 
the telephotometer expression is essentially identical to equation (10) . 
The differences between the two forms (eqs. (10) and (B5)) are confined 
to numerical constants affecting the magnitudes of the integral scale and 
intensity of fluctuations deduced from the data. The interpretations 
assigned to the x and y axis intercepts are the same for both cases 
(Markoff ian and Gaussian) . For a given set of data, such as any one of 
the configurations shown in figure 5, the Gaussian interpretation of the 
intercepts leads to values that differ by factors of 2 (approximately) 
from the corresponding values from the Markoff i an assumption: 


l 


Gauss 


2l Mark 



An experiment designed to measure the integral scale and intensity simul- 
taneously by an independent method is required to show which assumption is 
the more applicable to the turbulent boundary layer. 


CONCLUDING REMARKS 


To summarize, it has been shown analytic all y that plane light waves 
traversing a turbulent boundary layer undergo forward scattering which 
introduces a reduction in the transmittance ratio. Further, it has been 
demonstrated that photometric data obtained by measurements taken in the 
scattered field are in functional agreement with the predictions of the 
analysis. The interpretation of experimental data in terms of the turbu- 
lence parameters, however, depends upon the turbulence model employed and 
is in doubt by at least a factor of 2 for the models studied. 

Thp> principal advantage of the optical method described here is that 
the phenomenon of scattering is dependent only on the density fluctuations. 
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The measurements, therefore, are not hampered "by a mixed response to 
several fluctuating variables as is the hot-wire anemometer. Other advan- 
tageous features are that no probe and a minimum of electronic circuitry 
is required. 

The greatest limitation of the method is that, inherently, the results 
obtained are integrated averages over the volume of turbulence irradiated 
by the incident beam. Also, the possibility exists that for conditions of 
high losses, the results can become severely contaminated by the effects 
of secondary scattering, which were ignored throughout the analysis. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif. , 


for Aeronautics 
Apr. 9, 1956 
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APPENDIX A 

RELATIONSHIP BETWEEN DIEIECTRIC FLUCTUATIONS 
AND DENSITY FLUCTUATIONS FOR AIR 


For the system of units force, length, time, and electromagnetic unit 
of charge (emu), the wave velocity in any. medium is given by (ref. 8) 


1 

«/Te 


= v 


(Al) 


The velocity of wave propagation in a vacuum is the speed of light 


1 

•Tv- 


c 


(A2) 


The refractive index of a medium is defined by the ratio of the speed of 
light in a vacuum to the speed of light in the medium: 

For gases, the permeability to electromagnetic radiation is essentially 
that for a vacuum; the dielectric constant, however, differs by small, 
measurable amounts. For gases, then, the refractive index is given to 
good approximation by: 


i\ = 



(A4) 


Squaring the logarithmic differential of equation (Ah) yields the relation 
between fluctuations in the refractive index and the dielectric constant 
of the medium: 



(A5) 


For air, the refractive index is given by 
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where C is the Gladstone-Dale constant, 0.117 ft 3 /slug at 5200 A. 
Comb i tv? ng equations (A5) and (a 6) gives the relation between dielectric 
fluctuations and density fluctuations for air: 



= ^C 2 (dp ) 2 


(A7) 
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APPENDIX B 

TELE PHOTOMETER EQUATION FOR THE CASE OP A 
GAUSSIAN CORRELATION COEFFICIENT 


Thp scattering cross section derived under the assumption of a 
Gaussian form of the correlation coefficient is (Wheelon and Muehmore , 
p. 1452, ref. 3) 


where 



(Bl) 


(B2) 


(The relation between the solid angle and the scattering angle is the 
identity given previously (eq. (5)) w/4jt = sin 2 (0/2)). For cases of for- 
ward scattering wherein essentially all of the scattered energy can he 
found within very nma.1 1 solid angles, the upper limit in equation (8) can 


he taken as 4it as before. 

(Bl) and (B2) and noting that 


using equations 
obtains 


Carrying out the integration of equation (8) 

e"^ « 


-1 


, % 


43 1 


one 


(B3) 


2k 


W 


<2y 


For the small 3 angles encountered in forward scattering measurements, the 
first two terms of the exponential expansion can be used to obtain: 



(B4) 


or 
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■sTit 2rt 2 Z 
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1 



+ 


1 K 



tan 2 


% 

2 


(B5) 


where the sine has "been replaced by the tangent to correspond to equa- 
tion (10) . Comparison of equations (B5) and (10) shows that the forms 
are the same; the difference between the two results lies in the constant 

coefficients -/jt/2 and 2\f~n that are found in equation (B5) . The term 
K is essentially equivalent to 1 + K because K » 1. 
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Figure 1.- Experimental arrangement for viewing a light beam scattered by turbulent density 

fluctuations . 
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Figure 2.- Schematic view illustrating that only those' rays that are 
scattered through angles that are less than the half-angle of the 
focal-plane aperture are received hy the telephotameter. 
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Principal scattering cone containing one- 



Figure 3 • - Principal scattering cone for an incident plane wave for the 
scattering cross section given hy Booker and Gordon (angles greatly 
enlarged) . 
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Figure 5 .- Demonstration of the linear dependence obtained for telephotometer measurements in the 
scattered light field for fixed free-stream conditions and variable boundary- layer thickness. 
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